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Abstract 

To understand the effect of molecular environment on the electronic and magnetic properties of 
the single- molecule magnet (SMM) Mni2, we explore two possible means for adding extra electrons 
to molecule. We explore both substitution of Mn ions by Fe ions and the inclusion of neighboring 
electronic donors. For both possibilities we calculate, within density-functional theory the elec- 
tronic structure, the total ground-state spin and ordering, the magnetic anisotropy barrier, the 
transverse magnetic anisotropy parameter E which is responsible for some measured tunneling, 
and the tilting angle of the easy axis from the z axis. Our calculations show that the total spin 
increases with increasing number of extra electrons except for the case of Mn§Fe4 where the re- 
sulting ground state has a low spin. The calculated energy gaps between the unoccupied and the 
occupied orbitals exhibit no clear trend as a function of number of extra electrons. The calculated 
magnetic anisotropy barrier decreases with increasing number of extra electrons and can be di- 
rectly chased to a quenching of Jahn- Teller distortions at the sites where the additional electrons 
are localized. The values of E and the easy-axis tilting angles for the geometries with one- and 
two-extra electrons are significantly larger than those induced by solvent disorder. 



PACS numbers: 75.50.Xx, 71.15.Mb, 75.30.Gw, 75.30.Et 



I. INTRODUCTION 



There has been a great demand for high-density magnetic recording media that can store 
a bit of information on a nanoscale magnetic particle or molecule. Periodic arrays of high- 
spin molecules (each molecule having a volume of a few nm 3 ) could be potentially useful for 
this purpose. These high-spin molecules consist of several transition metal ions which are 
strongly coupled to each other and are surrounded by various ligands. The intermolecular 
interactions are quite weak and each molecule has a large energy barrier to magnetization 
reversal (i.e. magnetic anisotropy barrier). This leads to a behavior analogous to that of 
single-domain magnetic nanoparticles in an external magnetic field at low temperatures. 
Because of their behavior, the molecules are referred to as single- molecule magnets (SMMs). 
Although typical SMMs are paramagnetic, they exhibited magnetic hysteresis loops at low 
temperatures^* 2 - similar to that of ferromagnets at room temperatures. This is because 
their magnetic moments are frozen so that they relax on scales that are slow compared to 
the measurement time. The striking observation on SMMs was that they revealed quantum 
tunneling of magnetic moment 3 through the magnetic anisotropy barrier (MAB), which were 
confirmed by various experiments including temperature independent magnetic hysteresis 
loops 1,2 and temperature independent relaxation time of magnetic moment 4 . 

Among hundreds of synthesized SMMs, [Mn 12 12 (CH 3 COO) 16 (H 2 0)4] 
•2(CH 3 COOH)-4(H 2 0) (hereafter Mn 12 ), 5 has been the most extensively studied ex- 
perimentally and theoretically. The strong interest is due to Mni 2 being the first 
synthesized SMM and to the fact that to date it provides the largest MAB of 65 K.^^ The 
MAB is directly related to the temperature below which magnetic moments are frozen. A 
general spin Hamiltonian for a SMM is, to lowest order, 

Ho = DS 2 Z + E(S 2 X - S 2 y ) + gfi B B ■ S (1) 

where S z is the spin operator projected onto the z axis which is the magnetic easy axis, D 
and E are the uniaxial and transverse second-order magnetic anisotropy parameters caused 
by spin-orbit coupling. The last term is the Zeeman interaction. In 7f , the odd-order terms 
in S vanish due to time reversal symmetry at zero field. High-field, high-frequency electron 
paramagnetic resonance (EPR) experiments on the SMM Mni 2 provided D = —0.55 K&LIL 
Considering the S4 symmetry of a Mni 2 molecule (without solvent molecules), we know that 
E — 0. Quantum tunneling of magnetic moment is then allowed for pairs of states that 
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satisfy a certain selection rule: The two states must have ^-eigenvalues differing by integer 
multiples of 4. However, tunneling between states which violate this selection rule are clearly 
observed in magnetic tunneling measurements-*^ To resolve this puzzle, three different the- 
ories have been proposed . 9 ! 10 ' 11 ' 12 Chudnovsky and Garanin 9 - proposed any dislocations in 
single crystals could produce a broad continuous distribution in E. Detailed x-ray diffraction 
studies by Cornia et a/. 10 suggested that disorder caused by acetic-acid (CH3COOH) sol- 
vent molecules could break the S4 symmetry and provide nonzero locally varying E values. 
Magnetic quantum tunneling^ 3 - and high-field, high-frequency EPR experiments 1 ^ showed 
that their experimental data are not in quantitative accord with the Cornia et al.'s model^ -. 
Recent density-functional theory (DFT) calculations 1 ^* showed that for quantitative compar- 
ison with experiment, geometry relaxation must be considered. Pederson et al. suggested 
that coupling between the magnetic and vibrational degrees of freedom leads to the standard 
fourth-order form of the anisotropy Hamiltonian for isotopically pure samples.— However, 
they showed that for isotopically impure samples^ there are two different vibrational effects 
that cause symmetry and a coupling between M and M ± 2 states and that the symmetry 
breaking would be random assuming the isotopic substitution was. 

Very recent experiments on the SMM [Mn 12 12 (0 2 CCH 2 Br) 16 (H 2 0) 4 ]-8CH 2 Cl 2 (hereafter 
Mni 2 -BrAc) showed that there is a broad distribution in internal transverse fields (7.3° for 
a standard deviation).— For this case the molecules of crystallization CH2CI2 have higher 
symmetry and probably do not break the intrinsic S4 symmetry in contrast to the acetic acid 
solvent molecules in the standard Mni 2 crystal. It has been suggested that the distribution 
in internal transverse fields could arise from a distribution in the magnetic easy axes of 
individual molecules . 16117 But the origin of the distribution in the magnetic easy axes remains 
unknown. In this regard, we consider, within DFT, the possibility that extra electrons when 
accepted by Mn 12 molecules, will exhibit a disorder- or defect-induced behavior that is 
observed in the experimental samples. 

The SMM Mni 2 has four Mn 4+ ions at the corners of the inner cubane and eight Mn 3+ 
ions in the outer crown (Fig. [TJ. The magnetic moments of the eight Mn 3+ ions (S = 2) 
were observed to be antiparallel to those of the four Mn 4+ ions (S = 3/2), which leads to a 
total ground-state spin of S — 10.— The observed magnetic structure and the MAB for the 
ground state of the Mni 2 has been calculated in recent DFT calculations.— The exchange 
interactions between metal ions have also been calculated 20 within DFT. Both the MAB and 
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exchange interactions are in very good agreement with experiments.— Despite an increasing 
number of synthesized SMMs, it is not well understood what types of variations in the the 
molecular environment influence the magnetic structure, the intramolecular exchange inter- 
actions, and the MAB of the molecule. Hendrickson's and Christou's groups synthesized 
two types of anionic SMM Mn 12 and performed magnetic measurements on them . 22 i 23 i 24 
One type is gleaned by substitution of four Mn ions by four Fe ions in the ordinary Mni2 
structure, [Mn8Fe40i2(CH 3 COO)i6(H 2 0)4] (hereafter experimental Mn 8 Fe 4 ).— The other is 
synthesized by adding one or two PPh 4 molecules (Ph=C 6 H 5 ) to a Mn 12 molecule contain- 
ing different ligands from the ordinary Mn 12 , (PPh 4 ) [Mni 2 Oi 2 (0 2 CCH 2 CH 3 )i 6 (H 2 0) 4 ] (here- 
after (PPh 4 )Mn 12 ) and (PPh 4 ) 2 [Mn 12 12 (0 2 CCHCl 2 ) 16 (H 2 0) 4 ] (hereafter (PPh 4 ) 2 Mn 12 ). A 
PPh 4 molecule has one unpaired electron with Aj symmetry and may be thought of as a 
molecular alkali . 23 > 24 Magnetic measurements on both types revealed that their total spins, 
magnetic susceptibilities, and MABs significantly changed with the number of extra elec- 
trons but that the two types are partially but not entirely homologous . 22 i 23 i 24 As far as we 
know, no theoretical studies on these interesting results have been reported. 

As implied earlier, our motivation is twofold. One is to understand how the microenvi- 
ronment within a Mn 12 molecule controls its magnetic properties. The other is to view extra 
electrons as disorder or defects in the Mni 2 structure and to understand their effects on S 4 
symmetry-breaking terms in the spin Hamiltonian Eq. ((TJ. Motivated by the experiments 
described in the previous paragraph, we employ DFT to investigate two categories of envi- 
ronments where extra electrons can be transferred to metal ions in a single Mni 2 molecule: 
(i) substitution of Fe for Mn and (ii) addition of electronic donors such as K to the molecule. 
In the first case, since an Fe atom has one more electron than a Mn atom, replacing Mn by 
Fe will play the role of donating one extra electron to the ordinary Mni 2 structure. In the 
second case, since a PPh 4 molecule has one unpaired electron, addition of PPh 4 molecules to 
a Mn 12 molecule is equivalent to adding single-electron donors to the Mn 12 molecule. In our 
calculations, K atoms were used instead of large PPh 4 molecules because a K atom also has 
one unpaired electron and it reduces computational cost. In the sense that there are extra 
electrons in the ordinary Mn 12 structure, we call the SMMs under the new environment the 
anionic SMM Mn 12 . For both categories, we consider one-, two-, and four- extra electrons 
donated to the Mni 2 molecule. For each anionic state, we calculate the electronic struc- 
ture, the total magnetic moment, and the magnetic anisotropy parameters, using DFT. We 
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find that Fe substitution provides similar results (not exactly same) to K addition and that 
some microenvironments drastically alter the magnetic structure. In Sec. II, we describe our 
detailed model and method. In Sees. Ill and IV, we separately present DFT calculations 
on the electronic and magnetic properties for Fe substituted and K added geometries. In 
Sees. V and VI, we present our discussion and conclusion. 

II. MODEL AND METHOD 

Our DFT calculations 2 ^ are performed with spin-polarized all-electron Gaussian-orbital- 
based Naval Research Laboratory Molecular Orbital Library (NRLMOL) 26 . Here the 
Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA) is used for the 
exchange-correlation potential. 27 For all atoms full basis sets are used with a fine numerical 
integration mesh unless otherwise specified.— Prior to geometry relaxation, an initial spin 
configuration must be assigned with an initial total moment that can either be fixed or 
optimized. An optimized geometry is obtained when forces on all atoms become small. To 
find the ground-state spin configuration, different spin configurations with different total 
magnetic moments are considered. The second-order magnetic anisotropy parameters are 
calculated for an optimized geometry using the technique described in Ref. jl9|. 

As initial geometries for the two types of environmental changes, we start with a sim- 
plified form of a Mni 2 molecule with S4 symmetry, [Mni 2 Oi2(HCOO)i6(H 2 0)4]. Here 16 
acetates (CH 3 COO) in the experimental Mn 12 were replaced by 16 formates (HCOO) (refer 
to Fig.Q). Earlier calculations on the Mn 12 molecule showed that the energy gap between the 
majority HOMO (highest occupied molecular orbital) and the minority (majority) LUMO 
(lowest unoccupied molecular orbital) is about 0.804 eV (0.438 eV).— The calculated ma- 
jority HOMO-LUMO energy gap agreed with the result of the recent electrical resistivity 
measurement, 0.37±0.05 eV£& The total ground-state moment was calculated to be 20/is, 
which agreed with experiment. The calculated second-order MAB was 54 K in zero field, 
which did not change much with our simplification of the molecule. Hereafter, unless spec- 
ified, Fe substitution and K addition will be carried out in the above simplified form of a 
Mni 2 molecule. 

For the optimized geometry of the Mn 12 molecule, the first four unoccupied majority or- 
bitals directly above the Fermi level are close to each other in energy but well separated from 
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the rest of unoccupied orbitals. The four orbitals contain the LUMO, an orbital 0.008 eV 
above the LUMO, and doubly degenerate orbitals 0.012 eV above the LUMO. The fifth or- 
bital is 0.212 eV above the LUMO. Because the third and fourth orbitals are degenerate, the 
case of three extra electrons is not considered in our DFT calculations. But our prediction 
of the case based on symmetry will be provided in Sec. V. We first examine whether there 
are any favorable sites for the extra electrons in the Mni2 geometry when four Fe ions are 
substituted in a S4 symmetric fashion. Then when four-, two-, and one-extra electrons are 
accepted in a Mn 12 molecule under each type of environment, we calculate HOMO-LUMO 
energy gaps, the total ground-state moment, and the ground-state spin configuration. In 
various ways, we examine whether the extra electrons are well localized at some of the Mn 
sites. Using the optimized geometry with the ground-state spin configuration, we calculate 
the MAB, the transverse magnetic anisotropy parameter E, and the magnetic easy axis. 
When experimental data are available, comparison between calculations and experiment are 
discussed. 

III. DFT CALCULATIONS: IRON SUBSTITUTION 

The S4 symmetry and the geometry of the Mn 12 provide three symmetrically inequivalent 
Mn sites such as one inner cubane site (a) and two inequivalent outer crown sites (J3 and 
5) shown in Fig. d To locate an energetically favorable site for accepting extra electrons, 
we consider the three S4 symmetric configurations in which four Mn ions are replaced by 
Fe ions at the three inequivalent sites, [Mn 8 Fe40i 2 (HCOO)i 6 (H 2 0)4] (hereafter Mn 8 Fe 4 ). 
When four Mn ions at the sites a are replaced, the configuration is called the geometry 
a. A similar definition is applied to the geometries (3 and 5. The geometry 5 is shown 
in Fig. d The energies of the three tetra-Fe substituted geometries are self-consistently 
calculated with the total magnetic moment fixed and with the spin configuration illustrated 
in Fig. |2£a) (only the direction of each Mn moment is assigned such as spin up or down). 
Then we calculate HOMO-LUMO energy gaps and check whether Fermi filling is satisfied 
for the three geometries with two different fixed moments (24/ig and 22/ig). If there exist 
appreciable HOMO-LUMO energy gaps and the calculated orbitals abide by the Fermi filling, 
then the system is regarded as being magnetically stable. As summarized in Table [IJ the 
energy is the lowest among the six possible configurations when four Fe ions occupy the site 
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5 given the total moment of 24/ie- Local moments around the Fe sites for the geometries f3 
and 5 with the total moment of 24/ie, are found to be 3.94/i# and 3.95/is within a sphere 
with 2.23 Bohr radius (Table Q}. These local moments are definitely larger than those at the 
other outer Mn sites, 3.59-3.64/^; for the Mn 8 Fe4 geometries or those at the outer Mn sites, 
3.55-3.62/is, for the Mn 12 geometry. But the local moment at the Fe site for the geometry 
a is greatly reduced to —1.55 \lb- This is because the initial spin configuration [Fig. 
does not permit the extra electrons to be localized in the cubane. The same trend appears 
for the local moments of the geometry with the total moment of 22/i#. Our finding of the S 
to be energetically favored is also evidenced in the site calculated density of states projected 
onto the three inequivalent Mn(3cf) orbitals for the Mn 12 geometry. As shown in Fig. |3J the 
preferred site 5 has the largest majority density of states directly above the Fermi level. This 
finding is also corroborated by x-ray crystallographic data on the experimental MngFe^— 
Henceforth we consider only cases where Fe is substituted for Mn at the site 5. 



A. Four Fe substitution: MngFe4 

Our initial study showed that the Mn§Fe4 geometry with the total moment of 24fiB 
has lower energy than that with the total moment of 22/ig. The measurement on the 
experimental Mn 8 Fe4, however, revealed that the static magnetic susceptibility for this SMM 
has a completely different temperature dependence from that of the Mni2 (S = 10), which 
led to the ground state of S = 2 for the former.— So following the method discussed in 



Refs. 
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33l we examine several low-spin states to search for the ground state. In this study, 
we consider collinear spin states only. For a particular total moment, there are numerous 
spin configurations. We thus select a few spin configurations to test, considering that the 
exchange interactions between the Mn ions and Fe ions through the oxygen ligands (i.e. 
superexchange interactions)^ prefer antiferromagnetic coupling. We calculate the energies 
for six selected spin configurations with M s = 2 and for five chosen configurations with 
M s = 4 where M s is an eigenvalue of S z . Some of the configurations and their calculated 
energies are specified in Table ITT1 and Fig. Efb). Our calculations show that the ground 
state may have the total moment of 8fiB {M s = 4) with the spin configuration indicated in 
the sixth row of Table HP and Fig. |2Jb). Apparently this does not agree with the original 
experiment 22 but a very recent experiment on [Mn 8 Fe40i2(02CCH2Cl)i6(H 2 0)4] suggested 
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that the ground state may have a spin of S = 4.— Notice that the CH2CI molecules in this 
compound take up the position of the CH3 molecules in the experimental Mn§Fe4 compound. 
In our calculations, we terminated the carboxyl groups with hydrogen atoms rather than 
CH 3 radicals in the Mn 8 Fe4 geometry. Other possible reasons, such as a highly correlated 
many-spin wavefunction, for this discrepancy will be discussed later in this section and in 
Sec. V. 

We discuss the electronic and magnetic properties for the lowest-energy structure of the 
Mn 8 Fe4 which has the total moment of 8/xb and the spin configuration specified in Table ITT1 
and Fig. Efb). The energy gap between the majority HOMO and the minority (majority) 
LUMO is about 0.664 eV (0.247 eV). So the ground state is magnetically stable. The local 
moments of the inner Mn (the site a), outer Mn (the site (3), and outer Fe ions (the site 5) 
are found to be 2.61 fis, 3.48 and -4.16 fiB, respectively, within a sphere of 2.23 Bohr 
radius. Since the wave functions of Fe and Mn ions overlap with those of other neighboring 
atoms, the sphere around each ion does not capture a full moment. When we compare 
the local moments with those for the Mni2 shown in the second column of Table HI we 
find that only the local moment of the outer Fe ion substantially increases. Based on the 
calculations, we speculate that the inner and outer Mn ions and outer Fe ions correspond to 
Mn 4+ (S = 3/2), Mn 3+ (S = 2), and Fe 3+ (S = 5/2), respectively. This gives the ground- 
state spin of 5 = 4x3/2 + 4x2-4x5/2 = 4. The total MAB is calculated to be 33 K, 
which is approximately 60% of the MAB for the ordinary Mni2, 54 K (Table ITTT|) . Since the 
ground state has S4 symmetry, there is no transverse magnetic anisotropy, E = 0. Using the 
method described in Refslllll^ we calculate the projected MABs and the local easy axes of 
the three different metal sites. As shown in Table HVl (the third column), the local magnetic 
anisotropy of the Mn ions at the sites (3 mainly contributes to the total MAB. In contrast 
to the Mn 12 geometry, the Fe ions at the sites 5 have little contributions to the total MAB. 
This is supported by the calculated local easy axes of the three sites. The local easy axis at 
the site a is nearly in the x-y plane. The site 5 has a local easy axis 65.8° away from the z 
axis (global easy axis), while the site (3 has a local easy axis only 8.5° away from the z axis. 
As illustrated at the bottom of Fig. Ql an Fe 3+ ion in a distorted octahedral environment 
has exactly half-occupied e g and t 2g 3d orbitals, which leads to an absence of Jahn- Teller 
distortion as would be the case for a Mn 4+ ion and a Mn 2+ ion. Thus the total MAB arises 
from single-ion anisotropy of the four Mn ions at the sites (3 only. 
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Since the Mn§Fe4 geometry has a very different ground state from the Mni2, it is worth- 
while to explore the exchange interactions between metal ions within a Mn 8 Fe4 molecule. 
For this purpose, we consider collinear spin configurations with M s = 5, 6, 7, and 8 in ad- 
dition to the previously discussed several configurations with M s = 2, M s = 4 and M s = 12 
fTable ITT|) . The technique to calculate the exchange interactions using DFT was described 
elsewhere . 20 ! 33 ^ 34 As illustrated in Fig. Efb), there are six different exchange coupling con- 
stants within a Mn 8 Fe4 molecule such as J\, J 2 , J3, </ 3 , J4, and J' 4 . Here J 3 (J 4 ) is slightly 
different from J3 (J^) because of a small difference in the bond lengths. To compare with 
other references, we take an average of the two slightly different exchange constants. The 
spin Hamiltonian considered for the calculation of the exchange constants is 



where Eq is the background energy. All spins are of Ising type and they are numbered as 
illustrated in Fig. I2fb). From the DFT-calculated energies of the fourteen configurations 
(Table ITTJ) . we extract the seven unknown parameters, E , Ji, J 2 , ^3, J3, J4, and J4 using a 
least-square- fit method. Table 13 exhibits our calculated exchange constants in comparison 
with experimentally estimated values and the exchange constants for the Mni2 geometry.^ In 
contrast to the Mn 12 geometry, the Mn 8 Fe4 geometry has strong antiferromagnetic interac- 
tions between the outer Mn and Fe ions (J4 between the sites f3 and 5), which are comparable 
to the interactions between the inner Mn and outer Mn ions (Ji between the sites a and 
(5 and J 2 between the sites a and 8). This tendency was also observed in the experimental 
estimates. For a quantitative comparison of the calculated values with experiment, more 
refined experiments which are sensitive to the exchange constants are required. With the 
calculated exchange constants, one can construct the isotropic Heisenberg Hamiltonian. The 
true ground state is obtained by exact diagonalization of the Heisenberg Hamiltonian . 20 i 33 i 34 
The size of the Hilbert space to diagonalize is 207 360 000. We do not further pursue the 
ground state of the system by exact diagonalization of the large Hamiltonian matrix. The 
MAB varies with spin configurations from 16 K to 33 K as shown in the last column of 
Table |H] This large variation was not found for the SMMs Mni 2 and Mn4 monomer . ^i 34 
We speculate that the reason may be related to the fact that for the Mn 8 Fe4 the considered 
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spin configurations have energies of the same order of single-electron excitations fTable lTTTj) . 
B. Double and single Fe substitutions: MnioFe2 and MnnFe 

There are two inequivalent possible ways to replace two Mn ions at the sites 5 by Fe 
ions. One way is to place two Fe ions in a twofold symmetric fashion ("trans"). The other 
is to place Fe ions in a unsymmetrized fashion ("cis"). With the spin configuration shown 
in Fig. Efa), we obtain an optimized total moment of 22/xg for the trans and cis geometries. 
The experiment on the SMM (PPh 4 ) 2 Mn 12 suggested that the ground-state may have a 
total spin of S = 10M So we consider several spin configurations with M s = 10 (trans) 
that might have lower energies than the energy of the configuration with M s = 11. A good 
candidate is the configuration obtained by flipping only two spins coupled by J± from the 
configuration with M s = 11 (Refer to Figs. and (b)) We find that the configuration 
with M s = 10 has an energy 0.014 eV above the energy of the considered configuration with 
M s = 11 (Table IVTj). Thus, the ground state has a total moment of 22/xb. We calculate 
the electronic and magnetic properties for the ground state. For the trans geometry, the 
energy gap between the majority HOMO and the minority (majority) LUMO is 0.075 eV 
(0.432 eV). The value of E is 0.045 K and the MAB is 31 K (Table UTTJ) . The calculated 
projected anisotropy on the sites a, /3, and 5 implies that only the site (3 contributes to the 
total MAB (Table ITVj) . For the cis geometry, the energy gap between the majority HOMO 
and the minority (majority) LUMO is 0.030 eV (0.427 eV). The value of E is 0.020 K and 
the MAB is 33 K. The magnetic easy axis is tilted by 19° from the z axis fTable ITTTj) . 
Notice that the value of E for the cis geometry does not cancel out. This contrasts with a 
negligibly small E value calculated for the n = 2 cis configuration in a Mn 12 molecule with 
disordered solvent . 10 ^ 1 The MAB for the cis geometry agrees with that for the trans within 
numerical uncertainty. This barrier is significantly reduced compared to that for Mni 2 due 
to an absence of the Jahn- Teller distortion in the two Fe 3 + ions and noticeable values of E. 

When we replace one Mn ion at the site S by Fe, with the spin configuration shown in 
Fig. |2fa), we obtain an optimized total moment of 21yU#. The experiment on the SMM 
(PPli4)Mni 2 suggested that the ground state has a spin of S — 19/2.— We thus consider 
a few spin configurations with M s = 19/2 that might have lower energies than that for 
the considered configuration with M s = 21/2. A good candidate for this is obtained by 
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flipping only two spins coupled by J\ from the starting spin configuration. We find that the 
configuration with M s = 19/2 has higher energy than the M s = 21/2. So the ground state 
has a total moment of 21/zg. The ground state has the energy gap between the majority 
HOMO and the minority (majority) LUMO of 0.043 eV (0.424 eV). The total MAB is 41 K 
and E = 0.03 K. The magnetic easy axis is tilted by 10.7° from the z axis. 

IV. DFT CALCULATIONS: POTASSIUM ADDITION 

Adding electron donors such as K to the Mn 12 molecule, we consider the following two 
choices for the positions of the donors. One is to place the donors at the centers of mass of 
the acetic acid solvent molecules in the Mni2 (Fig. HJ. In this case, the donors are closer to 
the site (5 than the site 5 by 0.92 A. The other is to place the donors close to the favorable 
site 8 in the x-y plane (Fig. EJ). In our calculations, we use the former (latter) choice for 
four (two and one) K addition. If K atoms were too far away from the Mni2 molecule, then 
the unpaired electrons of the K atoms would not be transferred to the Mni2 molecule. The 
optimum distance between the K atoms and the Mn 12 molecule is found by relaxation of 
the geometry. When two K atoms are added, for direct comparison with experiment,— we 
consider the twofold symmetric Mn 12 molecule (Fig. EJ) as well as the S4 symmetric Mn 12 
molecule. The twofold symmetric Mn 12 geometry is obtained from the S4 symmetric Mn 12 
molecule by exchanging formates (HCOO) with neighboring water molecules in a twofold 
symmetric fashion (Fig. |SJ). All K added geometries we considered are charge neutral. We 
use a fine integration mesh and full basis sets for all atoms except for K atoms. We present 
results for four-, two-, and one-extra electrons added to the Mni 2 molecule. 

A. Four K addition: K 4 Mni 2 

Four K atoms are added to the Mn 12 molecule in a S4 symmetric fashion shown in Fig. 0]) . 
To compare with the experimental result, similarly to the four-Fe substituted geometry, 
we consider the three different spin configurations given in Table IVHI with three different 
magnetic moments (M s =2, 4, 12). The three considered configurations are speculated 
lowest-energy states for a particular M s . We find that the configuration with M s = 12 
provides the lowest energy among the three (Table EH]). In contrast to the Fe substitution, 
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when four K atoms are added, the high-spin state with S = 12 is the ground state. Plausible 
causes for this difference will be discussed later in this section and Sec. V. For the ground 
state, the local moments of the three sites (a, /3, and 5) captured within a sphere of 2.23 
Bohr radius are -2.56 3.79 and 4.16 respectively. Comparing these with the 
local moments for the Mn 12 geometry (Table IV11I vs Table HJ) confirms that the spin density 
of the extra electrons are transferred from the K atoms to primarily the site S in the Mni2 
molecule even though the site j3 is closer to the K cation than the site 5. This is in accord 
with the calculated projected density of Mn(3<i) orbitals (Fig. EJ). As specified in Table IVTTl 
the local moment of the site j3 for the M s = 12 state is 0.3 \xb higher than those for the 
M s = 2 and M s = 4 states. This indicates that for the M s = 12 state the spin density of 
the extra electrons is not entirely localized at the site 5 but there is some leakage to the site 
(3. This may contribute to the different ground state from the Mn§Fe4 geometry. 

We discuss the electronic and magnetic properties of the ground state. The energy gap for 
the majority HOMO and the minority (majority) LUMO is 1.04 eV (0.62 eV). See Table IvTTTI 
The calculated total MAB is 24 K. The projected MAB for the site 5 is calculated to be 
0.71 K and the easy axis of the Mn ion at the site is 51.1° tilted away from the z axis, which 
is qualitatively different from the projected MAB and the easy axis for the same site in the 
Mni2 geometry. See Table ITVl This verifies that the absence of a Jahn- Teller instability at 
the site 5 reduces the total MAB down to 24 K. This is only 44% of the total (2nd-order) 
MAB for the Mni2. We notice that the projected MAB of the site f3 is almost half of that 
of the same site for the Mni2 or Mn§Fe4 geometry. This might be attributed to incomplete 
localization of the spin density at the site 8. 

B. Double and single K additions: K2MJI12 and KM1112 

We first consider the case that two K atoms are located close to the site 5 where two water 
molecules are bonded in a twofold symmetric Mn 12 molecule. The geometry is referred to as 
K 2 Mn 12 (trans)*. The optimized total moment is found to be 22 //g. Similarly to the double- 
Fe substituted geometry, we find that the lowest-energy of the M s = 10 configurations is 
higher than that for the considered M s = 11 configuration as shown in Table IIXI So the 
ground state has a total spin of S — 11. The captured local moments on the Mn ions at the 
site S are found to be -4.15 fis within a sphere of 2.23 Bohr radius, while those of the rest 
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of Mn ions do not change much compared to those for the ordinary Mni 2 . The energy gap 
between the majority HOMO and the minority (majority) LUMO is 0.74 eV (0.31 eV). The 
calculated MAB is 31 K and E = 0.087 K (Table \VTH$ . 

Now we examine the case when two K atoms are added close to the sites 5 in the x-y 
plane to the S4 symmetric Mn 12 . As in the case of the Mn 10 Fe2 geometry, this molecule also 
has trans and cis geometries. For the cis, however, the HOMO-LUMO energy gaps do not 
open up upon geometry relaxation. Thus we consider the trans geometry only. The energy 
gap between the majority HOMO and the minority (majority) LUMO is 0.82 K (0.29 K). 
The calculated MAB is 34 K and E = 0.076 K. In comparison with the K 2 Mni 2 ( trans)*, the 
HOMO-LUMO energy gaps, the total MAB, and the E values agree with each other. 

When only one K is added close to the site 5 in the x-y plane, the optimized total moment 
is 21 fiB- Similarly to the single-Fe-substitution case, we find that the lowest-energy state 
with M s = 19/2 has higher energy than the considered configuration with M s = 21/2. So 
the ground state has a total spin of S = 21/2 (Table ITXjl . The ground state has an energy 
gap between the majority HOMO and the minority (majority) LUMO of 0.47 eV (0.036 eV). 
The calculated MAB is 46 K and E = 0.033 K. The magnetic easy axis is tilted by 4.1° from 
the z axis fTable IVIIIj) . 

V. DISCUSSION 

When several Fe ions were substituted for Mn ions at the most favorable 5 sites in the Mn 12 
geometry, we found that the total spin changes as a function of the number of substituted 
Fe ions (extra electrons) (See Table ITTTj) . The change of the total spin was the most drastic 
for the case of four Fe substitution. For the Mn 12 , MnnFe, and Mn 10 Fe2 geometries, the 
antiferromagnetic exchange interactions between the sites (3 and 5 (J4 and J' 4 ) are much 
smaller than J\ and J 2 [Fig. Efb)]. So the ground state has a high spin. Meanwhile, for 
the Mn 8 Fe 4 geometry, the values of J 4 and J' 4 are large enough to flip the spins of the 
inner Mn ions and the spins of the Fe ions in comparison to the spin configuration of the 
Mn 12 [Fig. |2fa)]. This leads to a ground state of a low spin. Although it is certain that 
the ground state must have a low spin, to determine the relative ordering of the S — 2 
and S = 4 states, we need to rely on exact diagonalization of the Heisenberg exchange 
Hamiltonian. The ground state may be a noncollinear quantum state (a linear combination 
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of Slater determinants). 

The MAB decreases with increasing number of extra electrons as shown in Fig. El because 
substituted Fe ions do not contribute to the total MAB due to a lack of Jahn- Teller distortion. 
Based on rigid band and intuitive argument, it was expected that the four Fe substitution 
would reduce the MAB more significantly than the MAB for two Fe substitution. This, 
however, was not found in our calculations. Although the reason is not entirely understood, 
it may be related to the result that the MAB varies up to a factor of 2 depending on the spin 
configurations given a particular geometry ( Table ITT]) . This tendency has not been observed 
for the Mni2 and the Mn4 monomer where different spin configurations gave almost the same 
MAB as that for the ground state . 20 ' 34 The major difference between the two cases is that for 
Mn 12 the spin flip energies are small compared to single-electron excitations, while for the 
Mn§Fe4 the spin flip energies are comparable or higher than the single-electron excitation 
energy. To further understand this large deviation, we calculate the projected MABs of the 
sites a, /3, and 7 for the Mn§Fe4 geometry with two different spin configurations (M s = 4 
and M s = 12). We found that the projected MAB of the site a for M s = 4 is smaller than 
that for M s = 12 and that the local easy axis of the same site for M s = 4 is not as close to 
the x-y plane as for the case of M s = 12 fTable llVj) . Since the projected MAB for the site a 
counteracts the total barrier or the projected MABs for the sites (3 and 5, a smaller barrier 
of this site would probably lead to a larger total barrier. 

To investigate the relationship 11 between the values of E for the single- and double-Fe 
substituted geometries (cis and trans), we use the spin Hamiltonian [Eq. (JTJ)] for the Mn i2 
and MnnFe and symmetry. Suppose that the spin Hamiltonians for the Mni 2 and MnnFe 
are 



where diagonalization of the matrices 7 ( -°* ) (7^) would provide the values of D and E. Then 
the spin Hamiltonian for the double-Fe substituted geometry is written, in terms of 7i[Mni 2 ] 
and W[MnnFe], as 



ft[Mn 12 ] 



E T^Vv, tt[Mn u Fe]= £ 7$^ 



V 1 



(3) 




ft[Mn 10 Fe 2 ] 



E ^[S,S V + R{S,)R{S V )\- e 7{2<W 



V 



(4) 




V ^(2t),(2c)o o 



(5) 



2 7 (i) T (2t) = ~(20 = 



(6) 
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(2c) = g (2c) = (1) _ (1) (2c) = (1) , (1) (J) 
Ixy ixz Ixz lyz > lyz Ixz 1 lyz V 1 / 

where R represents a proper symmetry operation and 7^ and y( 2c ' stand for the mag- 
netic anisotropy matrices for the trans and cis geometries, respectively. Here the magnetic 
anisotropy matrix 7 can be calculated within DFT.— Notice that a value of E may be de- 
termined by all of the off-diagonal elements of the matrix 7W, while the tilting angle of 
the magnetic easy axis is influenced only by their xz and yz components. For the single-Fe 
substituted geometry the value of E is obtained by the xy, yz, and xz components of the 
matrix 7W. The substantially larger tilting angle of the easy axis implies that the values of 
7W and 7W are not as small as those calculated for the Mn 12 with solvent disorder . 10 1 11 For 
the trans geometry, due to symmetry, 7^ = and 7^ = . Symmetry allows the nonzero 
off-diagonal element 7^ to depend on 7^ only [Eq. (j§J)]. So the value of E for the trans 
geometry is not twice as large as the value of E for the single-Fe substituted geometry. For 
the cis geometry, although the value of 7^ vanishes due to symmetry, there are nonzero 
values of y^f* and 7^ . Since y^f 1 and 7^ are not small, the easy-axis tilting and the value 
of E are substantial for the cis geometry ( Table UTT]) . The values of E calculated from the 
relaxed trans and cis geometries agree well with those obtained from Eqs. (jBJ) and (jZJ). The 
value of E for the cis is comparable to that for the trans or single-Fe substituted geometry, 
while the value of E for the n = 2 cis in the Mni2 with solvent disorder was an order of 
magnitude smaller that that for the n = 1 or n = 2 trans. 

When several K atoms were added to the Mn 12 geometry, we found that the magnetic 
structure and properties are primarily the same as those for Fe substitution except for a few 
points which will be described (Table IVIIII vs Table UTT]) . In both cases, the ground state 
has the same total spin for each number of extra electrons, except for four extra electrons. 
Unlike four Fe substitution, the four-K added geometry has the ground state of a high spin, 
S = 12. Similar to the case for Fe substitution, the MAB decreases with increasing number 
of added K atoms. Except for four-extra electrons, the K added geometries have the similar 
MABs to the Fe substituted geometries. The relationship between the values of E for the 
single- and double-Fe substituted geometries also holds for the single- and double-K added 
geometries. The values of E and the tilting angles of the easy axes are quite different for the 
two cases. We speculate that the difference between the two cases arises from incomplete 
localization of the spin density of the extra electrons on the sites 5 for K addition. Another 
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difference is that when one- or two-extra electrons are added, for K addition the energy gap 
between the majority HOMO and minority LUMO is an order of magnitude larger than that 
for Fe substitution. The calculated total MAB and the local easy axis for the sites S for 
both the single-K added and single-Fe substituted geometries, suggest that the geometries 
for one-extra electron have features in accord with the experimentally observed fast-relaxing 
species in a variety of low-symmetry SMMs Mni2 - 36l37l38l39l4 ° 

Although the case of the three-extra electrons was not considered in our DFT calculations, 
it is possible to calculate the values of D and E using the matrices 7^ and [Eq. Q] 
and symmetry. It is found that for triple-Fe substitution D = —0.20 K and E = 0.016 K, 
while for triple-K addition D = -0.29 K and E = 0.032 K fTables UTTl and IVTTTJ) . For 
completeness, the values of D and E for the K 2 Mn 12 (cis) is also calculated using Eq. (jSJ) 
and symmetry: D = —0.38 K and E = 0.003 K. The value of E for the K^Mn^cis) is an 
order of magnitude smaller than that for the MnioFe2(cis). In the event that such a system 
is experimentally realized, the calculated values of D and E can be compared to experiment. 



VI. CONCLUSION 



To examine the effect of extra electrons on the electronic and magnetic structure and 
the magnetic anisotropy, we considered, within DFT, two classes of molecules with extra 
electrons in the ordinary Mni2 geometry: Fe substitution for Mn and K addition as electron 
donors. Four Fe substitution implied that the sites 5 are energetically favorable for Fe 
ions, which agreed with experiment. Our calculations showed that the MAB decreases with 
increasing number of extra electrons, which is caused by absence of Jahn- Teller distortion at 
the sites where the extra electrons are localized. This tendency agreed with experiment and 
was confirmed by the calculated projected anisotropy on the sites. The total spin increases 
with number of extra electrons, except that the Mn§Fe4 geometry has the low-spin ground 
state due to the enhanced exchange interactions between the Fe and outer Mn ions. The 
S4 symmetry-breaking geometries such as the one- and two-extra electron geometries have 
E values that are at least by a factor of 3 larger than those associated solvent disorder. 
These larger E values will significantly expedite the tunneling and eventually reduce the 
total MAB. Additionally, the tilting angles of the easy axes are an order of magnitude larger 
than those for the Mni2 with solvent disorder, which will provide large internal transverse 
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fields that will further expedite tunneling. 
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TABLE I: Total fixed moments (2M S ), calculated local moments (sphere radii of 2.23 Bohr) asso- 
ciated with the three metal sites (a, (3, and 5 shown in Fig. [J), fip, and ^s, energy differences 
from the geometry 5 with M s = 12, AE, and whether the Fermi filling is satisfied, for the or- 
dinary Mni2 geometry and for the three initial tetra-Fe substituted geometries (a, (3, and 5), 
[MngFe40i2(HCOO)i6(H20)4], with the spin configuration of Fig. Ufa). In the geometries a, (3, 
and 5, four Fe ions are located at the four sites f3, 5, and a shown in Fig. ^ respectively. For 
example, Fig. ^ shows the geometry 5. 





Mn 12 


a 


P 


5 


a 


P 


8 




10 


12 


12 


12 


11 


11 


11 


Ha 


-2.58 


-1.55 


-2.59 


-2.54 


-1.84 


-2.59 


-2.57 




3.62 


3.60 


3.94 


3.59 


3.51 


3.74 


3.45 


HS 


3.55 


3.52 


3.64 


3.95 


3.48 


3.55 


3.83 


AE (eV) 




1.119 


1.186 





1.826 


1.106 


0.157 


Fermi filling 


Yes 


Yes 


No 


Yes 


No 


Yes 


No 
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TABLE II: Total moments (2M S ), down spins (the rest of spins are in an up state and each spin 
is numbered following Fig. 12(b)), Ising energy expressions [Eq. ©], the energy differences from 
the lowest-energy spin configuration (M s = 4) AE, the MABs for the fourteen collinear spin 
configurations in the tetra-Fe substituted geometry shown in Fig. ^ Here M s = 2 (V-2) denotes a 
spin configuration with M s = 2 which is different from the other three configurations with M s = 2. 
The exchange coupling constants Ji, J 2 , J 3 , J 3 , J 4 , and J 4 , are defined in Fig. 12(b). 



M s 


down spins 


Ising energy 


AE (eV) 


MAB (K) 


2 


2,4,5,6,7,8 


Eq - 9J 3 + 4.5 - 20 J 4 - 20J^ 


0.2764 


26 


2 (II) 


1,3,4,6,8,10 


Eq - 6J1 - 15 J 2 - 10 J 4 + 10 J' A 


0.4016 




2 (IV) 


1,8,9,10,11 


Eq - 15J 2 - 204 


0.2747 


29 


2 (V-2) 


5,6,7,11,12 


E Q - 6 Ji + 9 J 3 + 4.5J 3 ' - 10 J 4 - 104 


0.3941 


22 


4 


9,10,11,12 


Eq + 12Ji - 30J 2 + 9J 3 + 4.5J 3 ' - 2OJ4 - 20J^ 





33 


4(1) 


1,2,3,4,6,8 


Eq - 30 J 2 + 9 J 3 + 4.5 Jg' 


0.2808 


22 


4 (II-l) 


1,3,6,11,12 


Eq - 6 Ji - 9 J 3 + 4.5 - 10 J 4 - 10J^ 


0.4131 


21 


4 (II-2) 


2,3,6,11,12 


Eq + 6J1 - 15 J 2 - 4.5J^ - IOJ4 - 104 


0.3262 


27 


4 (III) 


1,6,7,8,12 


Eq - 12 Ji + 15J 2 - 2OJ4 


0.5265 


16 


5 


8,9,10,11 


Eq + 6 Ji - 15 J 2 + 9 J 3 + 4.5J^ - 204 


0.3312 


32 


6 


1,10,11,12 


Eq + 6 Ji - 15 J 2 - 10 J 4 - 10 J 4 


0.3407 


27 


7 


1,2,4,7,8 


£ - 6J1 - 15 J 2 


0.4240 


22 


8 


5,6,7,8 


Eq - 12Ji + 30J 2 + 9J 3 + 4.5J 3 / - 20J 4 - 20J^ 


0.5334 


24 


12 


1,2,3,4 


Eq - 12Ji - 30J 2 + 9J 3 + 4.5J3' + 20J 4 + 20^ 


0.6013 
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TABLE III: Calculated total spins, the second-order magnetic anisotropy parameters (D and E), 
easy-axis tilting angles 9 (angles between the local easy axes and the z axis), MABs, the energy 
gaps between majority HOMO and minority (majority) LUMO for the Mni2 and the Fe substi- 
tuted geometries. Except for MngFe3 the lowest-energy spin configuration has been used for each 
geometry. For MngFe3 results are determined from Eq. Q and the use of symmetry. 





Mnia 


MnnFe 


MnioFe2 (trans) 


MnioFe2 (cis) 


Mn 9 Fe 3 


Mn 8 Fe 4 


s 


10 


21/2 


11 


11 




4 


D(K) 


-0.54 


-0.42 


-0.29 


-0.33 


-0.20 


-2.10 


E(K) 


0.0 


0.030 


0.045 


0.020 


0.016 


0.0 


(deg) 


0.0 


10.7 


0.0 


19.0 




0.0 


MAB (K) 


54 


41 


31 


33 




33 


maj H - min L (eV) 


0.804 


0.043 


0.075 


0.030 




0.664 


maj H - maj L (eV) 


0.438 


0.424 


0.432 


0.427 




0.247 



TABLE IV: Calculated projected (local) MABs Aj and tilting angles of the local easy axes (pi for 
the sites, a, (3, and 5 (Fig. [TJ), in the Mn i2 (S = 10), Mn 8 Fe 4 (S = 4), Mn 8 Fe 4 (S = 12), Mni Fe 2 
(S = 11), and K 4 Mni2 (S = 12) geometries. 





Mn 12 (5 = 10) 


Mn 8 Fe 4 (S = 4) 


Mn 8 Fe 4 (S = 12) 


Mni Fe 2 (S = 11) 


K 4 Mni2 (S = 12; 


A Q (K) 


0.87 


0.21 


0.96 


1.16 


0.85 


A/3 (K) 


9.7 


9.5 


9.8 


10.4 


5.5 


As (K) 


9.3 


0.82 


0.41 


0.95 


0.71 


ip a (degree) 


83.8 


58.9 


87 


83 


83.4 


tpp (degree) 


9.8 


8.5 


8.3 


9.4 


9.1 


(ps (degree) 


36.2 


65.8 


61.5 


63 


51.1 
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TABLE V: DFT-determined and experimental exchange constants for the tetra-Fe substituted 
geometry shown in Fig. ^ and the Mni2 geometry. In this study a fine mesh and full basis sets were 
used for both the MngFs4 and Mni2 geometries. For the MngFe^ the fourteen spin configurations 
were used (Table [nj , while for the Mni2 , the eleven spin configurations were used (Ref. 2£jj ) . For 
comparison with other published results, a DFT-determined value of J3 ( J4) represents an average 
between the values of J3 and J3 (J4 and J4) shown in Fig. Efb). The exchange constants for the 
Mni2 shown here are slightly different from those reported in Ref. 20] because in the previous 
study smaller basis sets and fewer exchange constants were used (J3 = J3 and J4 = J4. 





Mn 8 Fe 4 (DFT) 


Mn 8 Fe 4 (Exp [Ref. 22]) 


Mni2 (DFT) 


Mnia (Exp [Ref. .21]) 


s 


4 


2 


10 


10 


Jl 


119 


432 


140 


119 


J2 


150 


397 


117 


118 


J 3 


7.2 


173 


7.2 


-8 


Ja 


122 


187 


24 


23 



TABLE VI: Total moments (2M S ), down spins, the energy differences from the M s = 21/2 (M s = 
11) geometry, and whether Fermi filling is satisfied, for the different spin configurations in the 
MnnFe (MnioFe2 trans) geometry. For the MnnFe geometry, an Fe ion occupies site 12 shown in 
Fig. EJb). For the MnioFe2 geometry, two Fe ions occupy sites 9 and 11 shown in Fig. Hb). The 
spin configurations for the M s = 21/2 and M s = 11 geometries are the same as Fig. Ha). The 
configurations for the M s = 19/2 and M s = 10 are obtained by exchanging spin 1 with spin 5 in 

Fig. lib). 





MnnFe 


MnioFe2 (trans) 




19/2 


21/2 


10 


11 


down spins 


2,3,4,5 


1,2,3,4 


2,3,4,5 


1,2,3,4 


AE (eV) 


0.035 





0.014 





Fermi filling 


Yes 


Yes 


Yes 


Yes 
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TABLE VII: Total moments (2M S ), down spins, local moments captured around the three sites, 
a, (3, and 5 (Fig. the energy differences from the lowest-energy spin configuration, whether 
Fermi filling is satisfied, and MABs, for the three different spin configurations in the four-K added 
Mni2 geometry, K4[Mni 2 Oi 2 (HCOO)i 6 (H20)4] shown in Fig.HJ For M s = 2, there are two slightly 
different local moments because of the twofold symmetry of the spin configuration. For M s = 4 
and M s = 12, the spin configurations are fourfold symmetric. 





2 


4 


12 


down spins 


2,4,5,6,7,8 


9,10,11,12 


1,2,3,4 


Ma (Mb) 


2.78, 2.52 


2.61 


-2.56 


M/3 (Mb) 


3.44, 3.48 


3.49 


3.79 




-4.20, -4.21 


-4.16 


4.16 


AE (eV) 


0.576 


0.359 





Fermi filling 


Yes 


Yes 


Yes 


MAB (K) 


20 


32 


24 



TABLE VIII: Calculated total spins, the second-order magnetic anisotropy parameters, easy-axis 
tilting angles, the MABs, the energy gaps between majority HOMO and minority (majority) LUMO 
for the K added Mni2 geometries. For K2Mni2 (cis) and KsMni2 results are determined from Eq. (|3"|) 
and the use of symmetry. 





KMn i2 


K 2 Mni 2 (trans) 


K 2 Mni2 (trans)* 


K^Mn^cis) 


K 3 Mni2 


K 4 Mni2 


s 


21/2 


11 


11 






12 


D(K) 


-0.46 


-0.36 


-0.35 


-0.38 


-0.29 


-0.16 


E(K) 


0.033 


0.076 


0.087 


0.003 


0.032 


0.0 


9 (degree) 


4.1 


0.0 


0.0 






0.0 


MAB (K) 


45.7 


34 


31 






24 


maj H - min L (eV) 


0.47 


0.82 


0.74 






1.04 


maj H - maj L (eV) 


0.036 


0.29 


0.31 






0.62 
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TABLE IX: Total moments, down spins, local moments, the energy differences from the lowest- 
energy spin configuration, and whether Fermi filling is satisfied, for two different spin configurations 
in the double K- and single-K added Mni2 geometries. In the trans geometry marked by *, the 
twofold symmetric Mni2 geometry (FigEJ) was used. 





K[Mn 12 0i 2 ] 


K 2 [Mni 2 0i 2 ]* {trans) 




19/2 


21/2 


10 


11 


down spins 


2,3,4,5 


1,2,3,4 


2,3,4,5 


1,2,3,4 


AE (eV) 


0.0823 





0.3225 





Fermi filling 


Yes 


Yes 


Yes 


Yes 
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FIG. 1: Schematic diagram of the fourfold symmetric tetra-Fe substituted Mni2 geometry, 
MngFe40i2(HCOO)i6(H20)4. The shading scheme of the atoms is as follows: darkness increases 
from H (white), to Mn, to Fe (largest spheres), to C, and to O (darkest). The three inequivalent 
metal ion sites are labeled as a, f3, and 5. The bottom figures illustrate occupancy of the e g and 
t 2g 3d orbitals for Fe 3+ (or Mn 2+ ) and Mn 4+ . 
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(a) (b) 




FIG. 2: (a) Initial spin configuration considered for the fourfold symmetric 
Mn8Fe40i2(HCOO)i6(H20)4 geometry to determine the most favorable sites for extra elec- 
trons. The inequivalent sites are indicated by the three different symbols, (b) Schematic diagram 
of the exchange constants, J±, J2, J3, J3, J4, and J4, where J3 and J3 (J4 and J4) are slightly 
different from each other because of a difference in the bond length by approximately O.lA. Each 
spin is labeled from 1 to 12. The sites a correspond to spins 1-4, the sites (3 to spins 5-8, and the 
sites 5 to spin 9-12. 
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FIG. 3: Calculated density of states projected onto the majority and minority Mn(3d) or- 
bitals at each of the three inequivalent sites (a, 0, and 5) for the fourfold symmetric 
[Mni20i2(HCOO)i6(H20)4] whose structure is exactly the same as Fig. fTJ with Fe replaced by 
Mn. The solid vertical line denotes the Fermi energy level. 



28 



FIG. 4: Schematic diagram of the fourfold symmetric K4[Mni20i2(HCOO)i6(H20)4] geometry, 
where two of the K cations are located slightly above the x-y plane and the rest two are slightly 
below the x-y plane. The shading scheme is the same as that for Fig. ^ except that the largest 
spheres are now K cations. The distance between the closest Mn (3 site (<5 site) and K is 4.403 
A (5.323 A). 
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FIG. 5: Schematic diagram of the twofold symmetric K2[Mni20i2(HCOO)i6(H20)4] geometry, 
where two K cations are close to the site 5 in the x-y plane. The shading scheme is the same as 
that for Fig. 0J The distance between the closest Mn and K is about 4.642 A. 
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FIG. 6: Calculated MAB vs number of extra electrons for the Fe substituted and K added Mni2 
geometries and the Mni2 geometry. 
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